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We study the multifractality (MF) of critical wave functions at boundaries and corners at the 
metal-insulator transition (MIT) for noninteracting electrons in the two-dimensional (2D) spin- 
orbit (symplectic) universality class. We find that the MF exponents near a boundary are different 
from those in the bulk. The exponents at a corner are found to be directly related to those at 
a straight boundary through a relation arising from conformal invariance. This provides direct 
numerical evidence for conformal invariance at the 2D spin-orbit MIT. The presence of boundaries 
modifies the MF of the whole sample even in the thermodynamic limit. 
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Anderson metal-insulator transitions (MITs), i.e., 
localization-delocalization transitions of noninteracting 
electrons, are continuous phase transitions driven by dis- 
order. At the transition, wave functions (WFs) are nei- 
ther localized nor simply extended, but are complicated 
scale invariant fractals exhibiting multifractal behavior 
characterized by a continuous set of scaling exponents 
[iIjSjBQI- In two dimensions (2D) the universal critical 
properties at a host of conventional phase transitions are 
known to be described by conformal field theories (CFTs) 

. It is natural to expect that disorder-averaged observ- 
ables at a localization transition in 2D are also governed 
by a CFT. If so, then conformal symmetry should im- 
pose severe constraints on averages of local quantities, 
including moments of WF amplitudes. 

In a recent Letter Subramaniam et al. extended the 
notion of multifractality (MF) to the boundaries of the 
sample ( "surface MF" ) , and showed that near boundaries 
critical WFs are characterized by MF exponents that are 
different from those in the bulk. Moreover, it was pre- 
dicted that the MF of the entire system depends cru- 
cially on the presence or absence of boundaries, even in 
the thermodynamic limit. In this Letter we study surface 
MF at the MIT in 2D for non-interacting electrons with 
spin-orbit scattering (symplectic universality class) 0, 
and extend the surface MF analysis to boundaries with 
corners ("corner MF" ). Conformal symmetry, if present, 
would lead (following [8|) to a simple exact prediction re- 
lating corner and surface MF exponents. Here, we show 
numerically that at the 2D spin-orbit MIT this predic- 
tion is indeed valid, thereby providing direct evidence for 
the presence of conformal symmetry at this MIT. We also 
confirm the dependence of the MF of the whole system 
on the presence of boundaries, as predicted in 

We begin by introducing corner and surface MF [^, 0| 
for a rhombus [Fig.[lja)] and a cylinder [Fig. [ijb)], both 
having edges of length L. All WFs ip{r) vanish at the 




FIG. 1: Systems studied: (a) A rhombus with the bulk, sur- 
face and corner regions of sizes I x I, I x h, and w x w sites, 
correspondingly; (b) A cylinder with the bulk {L x I) and 
surface (L x h) regions. 



boundaries. We define the corner (9) with opening angle 
6, surface (s), and bulk (b) regions of the rhombus, and 
similar regions s, b of the cylinder, as illustrated in Fig.[T] 
In each region 6, s, or b, the MF of WFs is characterized 
by the scaling of the moments of |-0(r)|^ with the system 
size L (all WFs ip{r) are normalized), 

L'i^J^Pi^ ^ L~^:^ ^ (x = 0,s,b,w), (1) 

where dx is the spatial dimension of each region (db = 
2, ds — 1, and dg — 0). The overbar represents the 
ensemble (disorder) average and the simultaneous spatial 
average over a region x surrounding the point r. r^, r^, 
and are the bulk, surface, and corner MF exponents, 
respectively. By x = w we label quantities computed by 
spatially averaging over the whole system (dw = 2) [6j]. 
Nonvanishing anomalous dimensions A^, 

A-^=T^-2q + d^, (2) 

distinguish a critical point from a simple metallic phase in 
which = 0. By the definition (Ullll), A^ vanish at q — 
and 1. The exponent /i defined in 6] is absent in Eq. 
^ because the local density of states is independent of 
energy at the spin-orbit MIT. The MF singularity spectra 
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f^ioi) are obtained from by Legendre transformation, 
= a^q - r,^ = dr^/dq. (3) 



f^{a) have the meaning of fractal dimensions: the num- 
ber of points rex, where scales as is pro- 
portional to L-f"'^°'\ This gives a direct relation between 
f^{a) and the distribution functions of WF amplitudes: 



-ln|i/;|Vlni. (4) 



Since f^{a^) ^ dg = 0, the ensemble average is essential 
[lo| for defining corner MF. 

Suppose that the q-th moment \ip{r)\'^i is represented 
by a local operator in an underlying critical field theory 
describing disorder averages. The scaling dimension of 
this operator will then equal Aq [11]. If the field theory 
possesses conformal invariance and if the operator is (Vi- 
rasoro) primary, then the relation = f between 
the surface and corner exponents can be derived 0] from 
the conformal mapping w = z^^^ . This yields 

a^-2=^(a^~2), = ^ [^(a^) - l] . (5) 

The validity of these relations provides direct evidence for 
conformal invariance at a 2D localization critical point, 
and for the primary nature of this operator. We note, 
however, that Eqs. ^ are valid only if a? > 0, because a 
is non-negative for normalized WFs [3, S] • K is expected 
[lot that for q > qe (where qg is a solution to = 0) 
the exponents become independent of q, while = 
[Eq. With the definition ^ this leads to a modified 
relation between and A^: 



(6) 



In [31 it was argued that when the MF in the whole 
sample with a smooth boundary is analyzed, the lowest 
of the Tq exponents for bulk and boundary dominates: 

— min(T^,r|). Points where the curves and 
intersect translate into linear segments on the plot of 
/^(a) (necessarily convex), interpolating between f^{a) 
and f^ia), see Fig. 3 in Q. 

In a recent paper Mirlin et al. pointed out that A^ 
obey the relation A ^ = A^'.^, which is expected to hold 
also for surface MF [12]. In two dimensions this leads to 



r(<) 



= 4-a^, (7) 



and implies that cannot exceed 4. 

To test all these theoretical predictions for the 2D spin- 
orbit MIT, we employ the "SU(2) model" defined in [H], 
which is a tight-binding model on a 2D square lattice with 
on-site disorder and a random SU(2) nearest-neighbor 
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FIG. 2: PDFs of logarithm of WF amplitudes on tori (red), 
in the bulk region of cylinders (black), in the surface region of 
cylinders (blue), and in the corner region with 6 = ti/2 (pink) 
and 6 = 7r/4 (light blue); L — 120. Inset: Semi- logarithmic 
plot. 



hopping. We consider four different lattice geometries: 
(i) torus, i.e., a square lattice with periodic boundary 
conditions (PBC) imposed in the x and y directions, (ii) 
cylinder [Fig. 1(b)] with PBC imposed in the x direction 
and open boundary conditions (OBC) in the y direction, 
(iii) square with OBC in the x and y directions, and (iv) 
rhombus [Fig. 1(a)] with 6 = tt/A and 9' = 37r/4 and 
OBC imposed in the x and y directions. In all these 
geometries the number of lattice sites is L^. 

For the scaling analysis the system size L is varied 
through L = 24, 30, 36, . . . , 120. For a fixed on-site disor- 
der strength Wc, we examined 6 x lO'* samples with differ- 
ent disorder configurations for each L. We have used the 
forced oscillator method ij] to diagonalize the Hamil- 
tonian, and extracted one critical WF from each sample 
which had the energy eigenvalue closest to the critical 
energy i?c = 1 at Wc = 5.952 (in the unit of hopping 
strength). For the results presented below, we have set 
I — L/6, /i = 1 for the cylinders [Fig. [IJb)], and w = 4 
for the corners [Fig. [Ua)]. We have numerically con- 
firmed that the exponents computed in the bulk regions 
of rhombi and cylinders agree with those of tori within 
statistical error bars. Also, the MF exponents for the 
surface region of rhombi are, within error bars, equal to 
those computed for the surface region of cylinders. In the 
following figures the bulk (surface) exponents are those 
computed for tori (for the surface region of cylinders). 

Figure [2] shows the probability distribution functions 
(PDFs) of In [■(/'(r)!^ measured for r at corners with angle 
9 — IT /A (light blue) and 9 — 7t/2 (pink), at the boundary 
of cylinders (blue), and in the bulk region of cylinders 
(black) at the fixed L = 120. Each PDF is normalized in 
the region where it is defined. The PDF calculated for 
tori is also shown in red, and it agrees quantitatively with 
the bulk PDF (black), as expected. Clearly, the PDFs 
for bulk, surface, and corner with 9 — 7r/2 and 9 = 7r/4, 
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FIG. 3: (a) Bulk (red), surface (blue), corner 6 — 71/2 (pink), 
and whole cylinder (orange) f{a) spectra, with error bars 
shown at integer values of q. Red, blue, and pink curves rep- 
resent f^(4 — a) + a ~ 2. Inset: Scaling plot of Eq. (O at 
q = 1 (filled circles) and q — 3 (open circles), (b) Corner 
/(a) spectra at 9 — 7r/4 (light blue), 7r/2 (pink), and 37r/4 
(green), with error bars shown at integer q (and at g = —0.5 
for 6 = 7r/4) . Curves represent the conformal relation ([5} . In- 
set: Numerical results for compared with Eqs. ([5]) (colored 
curves) . 



are all different, and, in this order, the peak position 
is shifted to the left, in agreement with the expectation 
that WF amplitudes should be smaller near edges. In 
the same order, the distributions become broader with 
longer (presumably power-law) tails at \4>\'^L'^ ^ 1. This 
means that for large q the moments \Tp[r)\'^i can become 
larger near edges (corners) than in the bulk, as the higher 
moments are dominated by long tails [ld| . 

We numerically obtain and f^{a^) from [see ([Tl[3|)] 



((InlV-P)). 



_ |V'(r)|29ln|?/'(r)|2 



-<lnL, (8) 



|i/;(r)|29 

^ [r (a^) - a^g - In L. (9) 

The inset of Fig.[3l^a) shows ((In I'i/'P})? as functions of L, 
computed for tori (x = b) , at the boundary of cylinders 
(x = s), and at the corners {9 = 7r/2) of squares. This 



inset exhibits distinct scaling behavior for bulk, surface, 
and corner regions for the displayed values of q = 1, 3. 

Figure[3ja) shows f^{a) of the bulk, surface, and cor- 
ner {9 — tt/2) regions. Clearly, in this order, the spec- 
tra f^{a) are seen to become broader and their max- 
ima are shifted to the right {a^ = 2.173 ± 0.001, 
al = 2.417 ± 0.002, a^/^ = 2.837 ± 0.003), in accor- 
dance with Fig. 2 and Eq. (|4]). (Recall that the maxi- 
mal values of f^{a) are the spatial dimensions d^-) The 
plot of /*(a) for the whole cylinder (Fig. [3l^a), orange) 
clearly represents the convex hull of f^{a) and f^{ct) 
[l,[l5[. Notice that f^{a) deviates from f^{a) already 
at f{a) « 1.5 (for q < 0). This confirms the prediction 
of (6( that the presence of boundaries drastically affects 
the MF of the system even in the thermodynamic limit 
and also in a typical sample (where /(a) > (To|). 

The data points of the bulk spectrum f^{a) [red dots 
in Fig. [^a)] lie on top of the red curve representing data 
points [l6| for f^{a\_q) in Eq. ^. This confirms Eq. ^ 
for the bulk which is the consequence of the symmetry re- 
lation I4I . Incidentally, the value of the typical bulk ex- 
ponent ag agrees with earlier calculations [l3,[I3| but not 
with [1^. Furthermore, our aQ satisfies 7r(ao — 2) = 1/Ac 
3, where Ac = 1.843 is the quasi-lD localization length 
at the MIT, normalized by the wire width, as obtained 
in pjj]. The surface spectrum P{a) (blue) is also seen to 
satisfy the relation ([7]) for 1 < a*' < 3, but there are dis- 
crepancies between the blue dots and the curve f^{oi\_g) 
when > 3 (g < -0.7) and < 1 {q > 2). Moreover 
it appears that can exceed 4, in contrast to < 4. 
This may question the validity of the symmetry relation 
of fl^ for boundaries, but we feel that computations on 
even larger system sizes and numbers of samples are nec- 
essary for drawing a definitive conclusion. 

Figure[3jb) shows the corner spectra f^{a) at = 37r/4 
(green), 7r/2 (pink), and 7r/4 (hght blue). As 9 de- 
creases, the peak position moves to the right {01^^'^ = 
2.558 ± 0.003, aj^ = 2.837 ± 0.003, and aj^ = 3.689 ± 
0.006) and the spectra become broader, indicating that at 
smaller the typical value of a WF amplitude is smaller 
but its distribution is broader. The numerical data (dots) 
are compared with the curves predicted from conformal 
invariance, Eq. using P{a^) of Fig. [3l^a) within the 
range 1 < ^3, where \q\ is sufficiently small to en- 
sure good numerical accuracy. The agreement between 
the numerical data and the predicted curves is excellent, 
confirming the presence of conformal symmetry. 

The inset of Fig. shows where the curves repre- 
sent computed with as input in Eq. ([5]). Note that 
a'' = 2 at g = 1/2 as a consequence of Eq. ([7|). We see 
that the numerical data for deviate from the predicted 
curves, Eq. (O, when < 1, in order to satisfy the con- 
straint > 0. We expect that in the limit L ^ 00, 
be given by Eq. ([5|) for q < qg and by = for q > qg. 

We note that the numerical results for exceed 4 
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FIG. 4: (a) The numerical data for the bulk exponents 
(red) are compared with their mirror image A\_q (grey). The 
green line represents ajj ~ 2. (b) The exponents Aq/[g(l — 
q)] for the bulk (red), the surface (blue), and corners with 
9 = 7v/4 (light blue), 7v/2 (pink), and 37r/4 (green). The 
curves represent the theoretical prediction, Eq. (|6}. Inset: 
Bulk (red), surface (blue), and corner (6^ — tt/4, light blue) 
exponents r^, and for a whole rhombus with 9 = n/A 
(black). 

when q < -0.1 for 9 = 7r/4, and q < -0.7 for 6 = tt/2 
[Fig-ISIb) and inset]. Even the maximum ag of f^{a) will 
exceed 4 for sufficiently small angles 0. On one hand, the 
maximum corresponds to q — Q where the numerics are 
most accurate. On the other hand, the maximal value 
f^{a) = dg = has a direct physical meaning as the 
dimension of the WF support, and must therefore appear 
on the f^{a) curve. Thus, our data strongly indicate that 
the symmetry relation of is violated for corners. 

The anomalous dimensions are computed numer- 
ically from |t/'(r)p'3/( |-0(r)|2 )'' ^ L^^'', which follows 
from Eqs. ^ and ([2|). FigurelD^a) shows the bulk anoma- 
lous dimension (red) and its mirror image across the 
q = 1/2 line A^*.^ (grey), both rescaled by q{l — q). 
Note that this rescaling magnifies small numerical errors 
around q — and q = I. Nevertheless the numerical data 
satisfy the relation A^ = A^'.^ of [H for -1 < g < 2 
where statistical errors are small. It is also clear from 
Fig. m^a) that A]^/[q{l — q)] varies with q, which means 
that the bulk spectrum f^{a) is not exactly parabolic. 

FigurelUb) compares A^ for bulk, surface, and corners 
with 9 = 7r/4,7r/2, and 37r/4. The solid curves represent 
the theoretical prediction ^ from the conformal map- 
ping, where A^ is taken from Fig. |4Ub). For sufficiently 
small values of |g| the numerical results of A^ are in good 
quantitative agreement with the prediction It is pre- 



cisely for small \q\ that the numerical data are most accu- 
rate (20| . This provides direct evidence for the presence 
of conformal symmetry at the 2D spin-orbit MIT. 

The inset of Fig.[4{b) shows the exponents for bulk, 

surface, and corners. We see that t^^^ (light blue) is 
constant for q > q^r/i ~ 1 reflecting the exchange between 

top and bottom lines in Eq. ^ which happens at aq^^ = 
0. It appears that r^^'' becomes smaller than both and 

for q > 2.5, which is when the corner exponent Tq^'^ 
controls the MF of the whole sample with a 7r/4 corner, 
as shown in black. In a sample without corners such as a 
cylinder, the surface exponent controls the MF of the 
entire sample for sufficiently large q. This confirms and 
generalizes the predictions made in Q. 

In summary, we studied bulk, surface, and corner mul- 
tifractality at the MIT in the 2D spin-orbit symmetry 
class. We provided direct numerical evidence for the pres- 
ence of conformal symmetry at this critical point, and 
confirmed predictions of Q that boundaries affect MF of 
the whole system even in the thermodynamic limit. We 
also tested the validity of the symmetry relation of 
for the bulk, surface and corners. It appears that the 
relation holds in the bulk, but is violated at corners. 
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